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Abstract. We prove that the critical temperature for the BCS gap equation 
is given by 

in the low density limit fj, ^ 0, with 7 denoting Euler's constant. The formula 
holds for a suitable class of interaction potentials with negative scattering 
length a in the absence of bound states. 



1. Introduction and Main Result 
The BCS gap equation [1] 



with E{p) — ^ [jp- — + |A(p)P, has played a prominent role in physics in the 
fifty years since its introduction. In ([1]), the function A is interpreted as the order 
parameter describing paired fermions (Cooper pairs) interacting via a local pair 
potential 2V(x), with Fourier transform V(p) = {^■k)~'^I'^ ^^^V (x)e~^'P'^ dx. The 
(positive) parameters T and /i are the temperature and the chemical potential, 
respectively. Originally, Eq. ([1]) was introduced to describe electrons in a crystal, 
where the interaction is mediated by phonons and is non-local. Special interest in 
local potentials of the form ([T]) has recently emerged for the description of ultra-cold 
fermionic gases (see, e.g., [2] and references therein). 

It was shown in [5] that the critical temperature for the existence of non-trivial 
solutions of the BCS gap equation can be characterized as follows. 

DEFINITION 1. Let > and G L^/^ be real- valued. Let Kt,^, denote the 
multiplication operator in momentum space 

- Ml 



^^^^ tanh(|p2 _ ^|/(2T)) ' 
The critical temperature for the BCS equation is given by 



= inf{r>0 : Xt.p + V^>0} (2) 
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It was proved in [5j that Eq. ([T]) has a non-trivial solution if and only if T < Tc. 
Note that this result gives a linear criterion for the existence of solutions of the 
nonlinear equation ([1]). An analysis of Tc for weak coupling (i.e., V replaced by AV^ 
with A ^ 1) was given in ^ (see also T;). In contrast, here we are interested 
in the low density limit fj, at fixed interaction potential V. In this regime, Tc 
turns out to be related to the scattering length of 2V, which can be conveniently 
defined as follows. 

DEFINITION 2. Let V e L^{R^) n L^/^{M.^) be real-valued, and let V{xy/^ = 
sgn(V^(a;))|V^(a;)|^/^. If —1 is not in the spectrum of the Birman-Schwinger operator 
then the scattering length of 2V is given by 




Note that since V £ L'^/^{W^) by assumption, the Birman-Schwinger operator is 
of Hilbert-Schmidt class by the Hardy-Littlewood-Sobolev inequality fS', Thm. 4.3]. 
Although it is not self-adjoint, its spectrum is real. The fact that —1 is not in its 
spectrum means that + V does not have a zero eigenvalue or resonance. Eq. ([3]) 
is the natural definition of the scattering length for integrable potentials. In the 
appendix, we shall explain why this definition coincides with the usual concept 
found in quantum mechanics textbooks. For an alternative definition allowing for 
more singular local behavior but assuming compact support of V and the absence 
of bound states see \TQ', Appendix A] . 

With the aid of Definitions [T] and [2] we can now state our main theorem. 

THEOREM 1. Assume that V{x){l + \x\) e L'^{R^) L^/'^{R^) is real-valued. 
Assume that the spectrum of V^^"^ ^\V\^^'^ is contained in (— l,c»), and that the 
scattering length a in is negative. Then the critical temperature 0^ satisfies 

lim (ln;^ + -^) =2-7-ln- (4) 

with 7 w 0.577 denoting Euler's constant. 
In other words, 

Tc^f^ (^^e^-' + o(l)^ e-/(2v^a) 

as /i 0. This formula is well-known in the physics literature [H [51 [TTJ \T2\ . 

The operator V^^^^\V\^^^ having spectrum in (—1, oo) implies, in particular, that 

p'^ + V does not have any bound states. The proof of Theorem [T] shows that the 
condition that y(a;)|a;| is integrable at infinity is indeed optimal; i.e., formula Q 
will in general not hold if V decays slower than Ixj"** at infinity. 

In the next section we shall prove Theorem [1] Our analysis of the operator 
KT,ti + 1^ in ([2| is similar in spirit to the spectral analysis of Schrodinger operators 
in two dimensions in [131 HI] ■ 

2. Proof 

Let — 1/A denote the smallest eigenvalue ofV-^^^^lV]^^^. By assumption, A > 1. 
It follows that for this value of A, -f XV > 0. 
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Since tanh< < min{l,<} for i > it is easy to see that Kt,^ > A ^(|p^ — 
/i| + 2T{\ — 1)) for any A > 1. With the choice of A as above, we conclude that 
Kt.^ + V > X-\p^ + XV - (1 + 2T(A - 1)) > A-i(-M + 2T(A - 1)). From this 
bound and ([2|) it follows immediately that Tc < /i/(2(A— 1)). Hence we restrict our 
attention to the case T < Dfi in the following, with D — 1/(2(A — 1)) > 0. Only 
the existence of such a D, and not its value, will be important in the following. 

According to the Birman-Schwinger principle (see [Sj Lemma 1]), Tc is deter- 
mined by the fact that for T ~Tc the smallest eigenvalue of 



i3^ = yi/2 1 |y|i/2 



equals —1. If the spectrum of Bt is contained in (— l,oo) for any T > 0, then 
Tc — 0. Alternatively, Tc is the largest T such that 1 + Bt has an eigenvalue 0. 
We decompose the Birman-Schwinger operator Bt as 

Bt = V^'^-^\V\"^ = V'''\\V\"' + m,{T)\V''^){\V\^'^\ + At^, , 



with 



(27r)3 7jj3 \Kt,^{p) p- 
Explicitly, At, p. is the operator with integral kernel 



We note that for small /i 

= ^ (in ^ + 7 - 2 + In ^ + o(l)) (5) 

uniformly in T for T < D^. This was shown in [6, Lemma 1]. 
Lemma 1. Under the assumption that |V"(x)|(l -I- |a;|)(ia; < oo, 

lim sup ^ .^J ^T.plb = 0- 

Oo<T<Dp fJ.'^'^m^[T) 

Here, || • ||2 denotes the Hilbert-Schmidt norm. 



Proof. Performing the angular integration, we can write the kernel of At,^ as 

27r^ Jo V P\x-y\ J \Kt,^{p) p^ J 

(6) 

We note that |6-isin6 - 1| < Cminj^^, 1} < Ch" for aU < a < 2 and 6 > 0. 
Hence we can bound, for any Z > 0, 



s\n.p\x — y\ ^ 



p\x - y\ 



< C 



p^Z^eiZ -\x- y\) + ^p\x~y\ei\x ~y\-Z) 0(2^ - 



+ C^p\x - y\e{p' - 2y) . 
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The kernel of At.^, is thus bounded by 

\AT,^{x,y)\ 



< |y(x)|l/2|y(y)|l/2^2 



c 



/2JI 



p dp 



+ \V{x)\'^'\Viy)\'^'\x-y\'^'ei\x-y\-Z) 



c 

1 



1 1 

KT,t,{p) P^ 

p'^/^dp 



p^/^dp 



It is not difficult to see that the p integral on the first line is bounded by (a constant 
times) fim^{T) for T < Dp,. Similarly, the integral on the second line is bounded 
by /i^/*m^(T). Finally, the last integral is bounded by p^^^. Since by assumption 
/g3 |y(a;)|(ix < oo and /^g 1^(2^)1 |a; — y\dxdy < oo, we conclude that 



lim sup — 



I 



0^ fi^^m^iT)' 



\A 



< c 



1/2 



\V{x)\\V{y)\\x-y\dxdy 



\x-y\>Z 



for some constant C > 0, uniformly in T for T < ZJ/i. Since Z was arbitrary, this 
proves the claim. □ 



Since 1 + V^^^p ^1^1^/^ is invertible by assumption, we can write 



Hence Tc is the largest T such that 



l + \/i/2p-2|y|i/2 



/2p-2|^|l/2 
Fl/2)(|^|1/2| + 



V^l/2)(|^|1/2| + 



A 



has an eigenvalue —1. According to Lemma [U the operator AT^^/m^ {T) is small 
for small /i, uniformly in T. Moreover, m^(T) is monotone decreasing in T and the 
rank one operator 



J _|^l/2\/|^|l/2| 

1 + ^1/2^-21^11/2 1" /Vl"! I 



has an eigenvalue Ana according to Definition [21 which is negative by assumption. 
Hence simple perturbation implies that 

lim m^{Tc) = ■ 
M-^o Ana 

In particular, T^, > for small enough fi. 

We conclude that m^{Tc) is order one as /i ^ 0. With the aid of Lemma [T] we 
thus see that for T — Tc, 1 + V^^^p~^\V\^^^ + ^t,^ is invertible for small enough 
II. This implies that Tc is determined by the fact that the smallest eigenvalue of 



1 + St = 



|yi/2)(|y|i/2| 
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is zero. Since the first factor does not have a zero eigenvalue, the second does for 
T = Tc- In other words, 

1 , 1 



^(|y|i/2| 



In order to evaluate this expression, we expand 
1 



(7) 



1 + V^/^^\V\^/^ + At, 



1 



1 



1 



" l + yi/2j^|y|i/2 i + yi/2^|y|i/2 ^'^l + yl/2j^|^|l/2 

1 1 1 

^ 1 + ^1/2 I V^l ^'''l + yl/2j^|y|l/2+A^^^ ^'^l + yl/2j^|^|l/2 • 

The first term, when inserted in ([71), yields l/(47ra). The last term term goes to 
zero faster than /i^/^ since, according to Lemma[Tl <C fJ,^''^. Inserting the 

second term into ^ yields (/|(sgnF)AT^^^|/) with 

l/)= . , ..uli..u. \V'^')- (8) 



l + 1 |y|l/2 



Lemma 2. Assume that f{x)\V{x)\'^^^{l + \x\) e L^{R^). Then 

1 



lim sup ■ , 



\{f\isgnV)AT,M)\^0. 



Proof. Using ^ we can write {f\{sgnV)AT,^i\f) as 



smp|x — y\ 



1 



1 



1 \ p'^dp 



27r2 



(9) 



dxdy . 



p\x-y\ 

Proceeding as in the proof of Lemma [l] but with a = 1 in place of a = 1/2 in the 
second term, we conclude that 



smp|a; — y\ 
p\x - y\ 



1 



< Z^C 



1 



1 1 

KtAp) 
1 



p dp 



X - y\e{\x - y\ - Z)C 
1 



p dp 
277 



1 



KtAp) p^ 



p^dp 



+ \x-y\'/^C 



V2^ 



1 



p^'^dp. 



The first and second integrals are bounded by ^im^iT) and y^TO^(T), respectively. 
The last integral is bounded by /x'^/^. Since Z is arbitrary, the same argument as 
in the proof of Lemma [T] leads to ([9]). □ 

If / defined in ([8]) satisfies the assumption of Lemma [2l we conclude that 
{f\{sgy\V)AT^,^\f) ^ /i^/^ as /i — > 0. In particular, this implies that 



lim ( 



47ra j 



Using IS]) this proves ([4]). 



6 



CHRISTIAN HAINZL AND ROBERT SEIRINGER 



Thus it remains to show that f{x)\V{x)\^/^{l + \x\) G L^{M.^), with / defined in 
. For this purpose, let g E i°°(R'^), and consider the expression 

1 



g{x)\V{x)\'/^l + \x\)f{x)dx = {g\V\'^'\{l + \x\) 



l + yi/2j^|y|i/2 



(10) 



We write 



1 + ^1/2 i|y|l/2 



1 - V 



1/2. 



1/2. 



For the contribution of 1 to (jTU]) we use 



{g\V\'/'\a + \x\)\V'/'] 



< 



\Vix)\{l + \x\)dx. 



(11) 



The remaining term is then bounded by 

\\g\V\'/'il + \x\m\2\\V'^% 



(1 



,1)1/2^1/2 1 



V 



1/2 



p. 



The first norm is bounded by ||5||oo(/ |V^(a;)|(l + The second is just 

1 /2 

. The first operator norm in the second fine is finite by assumption. Finally, 
the last operator norm can be bounded by the Hilbert-Schmidt norm, whose square 
is 



(i + i^i)in^)ii 



1 



F - y\ 



:\Viy)\dxdy<C\\Vh/2\\il + \x\)Vh/2 



The last inequality follows from the Hardy-Littlcwood-Sobolev inequality [9, Thm. 4.3] ; 
the right side is finite by our assumptions on V. This shows that (fTO|) is bounded 
by (a constant times) Uglloo, and hence f{x)\V{x)\^/'^{l + \x\) e L^. This completes 
the proof of Theorem [TJ 

Appendix A. Scattering Length 

In the following, we shall explain why the definition ^ of the scattering length 
coincides with the usual textbook definition. There, one considers the solution ip 
of the zero-energy scattering equation 

~Aip + Vip = 

subject to the boundary condition Ihn^^^^oo tplx) = 1. The scattering length a is 
then determined by the asymptotic behavior ipi^) ~ 1 ^ ct/l^^l for large For a 
rigorous definition for potentials with compact support and in the absence of bound 
states see [ini Appendix A]. 

Assuming that we have such a solution ip with ip — I ^ a/\x\ for large |x|, we 
have 

^Ai/' = ^A(V'-l) = l-V'- 



P 



p. 



This can be rewritten as 



or 



V 



l/l/2ly^ = yl/2(l_^)^ 
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Moreover, / Vip = / Atp — Ana, which foUows from integration by parts, using that 
I'^tpl « for large 

Under the assumption that 1 + V^^'^p^^\V\^/^ is invertible, we thus have 

Hence the rank one operator 

1 



1 + V^i/2j_|y|i/2 4;ra 



has an eigenvalue 1, which implies that 

/|1/|1/2| L |yl/2\ ^ 



1/1/2 J,|y|l/2 
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